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^ ! Abstract 

^ ■ 
Q^l The effect of tlie recently observed 2.6 a deviation of tlie muon anomalous magnetic 

p ■ moment (a^ = {g^ — 2)/2) from its Standard Model prediction is examined within 

^ ■ the framework of supergravity models with grand unification and R parity invari- 

• • ! ance. The constraints of the Higgs mass bounds, the 6 — *■ S7 bounds (including the 

.!^ I large tan/5 NLO corrections) and the cosmological relic density of light neutralinos 

r> I (including all slepton neutralino coannihilation effects) are included in the analy- 

c^ ' sis. For universal soft breaking, the Higgs and 6 — *> S7 bounds puts a lower bound 

mi/2 ~ 300 GeV, most of the parameter space now falling in the co-annihilation re- 
gion. The 2 a lower bound on the magnetic moment anomaly places an upper bound 
of mi/2 ~ 800 GeV. It is seen that mSUGRA requires that a^ ~ 50 x 10~^°. One 
finds for m^ > 114 GeV, that tan/3 > 5(7) for Aq = 0(-4mi/2) and for rrih > 120 
GeV, one has tan/3 > 15(10) for Aq = 0(— 4mi/2). The sparticle spectrum is now 
much constrained, and the reaches of the Tevatron RUN II, NLC, and LHC for new 
physics discovery are discussed. Dark matter detection rates are examined, and 
it is seen that future detectors now would be able to scan most of the parameter 
space. Models with non-universal soft breaking in the Higgs and third generation 
of squarks and sleptons are exmained , and it is seen that a new Z s-channel anni- 
hilation of neutralinos in the early universe is possible with dark matter detection 
rates accessible to the next round of detectors. 



1 Introduction 

Recently, the Brookhaven E821 experiment has measured the muon anomalous mag- 
netic moment, a^ = (g^ — 2)/2 for the yU^ to remarkable accuracy cD: 

a^ = 11659202(14)(6) x 10"^° (1) 

This had led to a 2.6 o" deviation from the prediction of the Standard Model (SM): 

a;^-P - af = 43(16) x 10"^° (2) 

Eq. (1) arises from an analysis of only about one quarter of the existing /i+ data. 
Further, the 2001 run for the fi~ should result in about three times as much data 
as has been analysed for the yU+ a) . Thus the statistical error in the measurement 
should decrease by about a factor of 2.5. In addition, the theoretical error in the 
hadronic vaccum polarization contribution should also be reduced (perhaps by a 
factor of 2) by the analysis of the Novosibirsk, DAPHNE and Beijing data. (For 
an analysis of the current theoretical accuracy of the hadronic vacuum polarization 
contribution see ^, and for an alternate view see o'.) Thus by the beginning 
of 2002 it should be possible to know whether the anomaly in the muon magnetic 
moment of Eq. (2) is real. If it is indeed real it would imply the discovery of new 
physics. Because of the importance of this possibility, we will here assume that the 
anomaly is real, and analyze some of its consequences. 

Supersymmetry offers a possible explanation of a deviation of a^ from its 
SM value, and initial calculations of a^ were made in the early 1980's using global 
supersymmetry a) . However, in unbroken global supersymmetry there is a theorem 
that the total value of a^ should vanish B*, i.e. 

aT + «f '^ = (3) 

Thus one needs broken supersymmetry to get a non-zero result, and how to do this 
in a phenomenologically acceptable way within the framework of global supersym- 
metry was problematical. The advent of supergravity (SUGRA), where acceptable 
spontaneous breaking of supersymmetry occurs, led to the SUGRA grand unified 
(GUT) models \^ . The first (partial) calculation of a^^^^"^ was done in l3J, and the 
first complete calculation in hJ. In the SUGRA GUT models, SUSY breaking at 
the GUT scale Mq, triggers the breaking of SU{2) x U{1) at the electroweak scale 
Mew, and hence the SUSY mass scale Msusy obeys 



'SUSY = ^WEW 



(H) (4) 



setting the scale of the SUSY masses to be in the 100 GeV - 1 TeV range, a range also 
needed to avoid the hierarchy problem. This mass scale was further supported by the 
fact that the LEP data is consistent with grand unification if the SUSY masses lie 
in the range ~ 100 GeV - 1 TeV. Further, SUGRA models with R parity invariance 
have a dark matter candidate, the lightest neutralino, Xi^ with the astronomically 
observed amount of relic density when the SUSY masses are in this range. Thus 
with the mass range implied by the SUGRA GUT models, it was possible to predict 
roughly the size of a^'^^^, and in fact in 1^'' it was argued that a^^'"'^'^ should 
become observable when the experiments had a sensitivity of about 20 x 10~^°, as 
is the case for Eq. (2). 

We consider here, then the effect of q^ugra ^^^^ SUGRA GUT models pos- 
sessing R parity invariance. In particular, we consider (1) the mSUGRA model 111) 
with universal soft breaking at Mq, and (2) models with non-universal scalar masses 
at Mg for the Higgs bosons and the third generation of squarks and sleptons. 

2 Calculational Details 

In the following we impose the current experimental bounds on the SUSY parameter 
space and list now the important ones: 

1. Accelerator bounds: (i) We consider two possible values for the light Higgs 
mass rrih > 114 GeV and rrih > 120 GeV. The first is the current LEP 
bound c^'' while the second represents a bound that the Tevatron Run IIB 
could achieve, if it does not discover the Higgs boson, (ii). We assume for the 
b —>■ sj branching ratio the roughly 2cr range: 

1.8 X 10"^ < BR{b -^ S7) < 4.5 x 10^^ (5) 

2. Relic density bounds for the lightest neutralino. We chose here the range 

0.025 < n^oh^ < 0.25 (6) 

where Q^o = p^o/pc, and pc is the critical density to close the universe: pc = 
SHq/SttGn (and Hq is the Hubble constant, Hq = h 100 km/sec Mpc, Gn the 
Newton constant). The lower bound in Eq. (6) is smaller than the conventional 
value of ~ 0.1, but can take into account the possibility that there is more 
than one species of dark matter. Our results here, however, are insensitive to 
this lower bound and would not change significantly if it were raised to 0.05 
or 0.1. 



3. We use a 2a bound on the muon magnetic moment anomaly of Eq. (2): 

11 X 10-1° < a^UGRA < 75 X ^^-x^ (7) 

As we will see these will combine to put strong constraints on the SUSY parameter 
space and allow us to greatly strengthen the predictions for the neutralino-proton 
dark matter cross sections o'j^o.p, as well as determine what parts of the SUSY mass 
spectrum can be seen at the Tevatron Run II, LHC and NLC accelerators. 

In order to get accurate results, it is necessary to include a number of 
corrections in the calculation, and we mention here the most important ones: (1) 
The allowed SUSY parameter space is now quite sensitive to the Higgs mass, and so 
an accurate calculation of m^ is important. We include here the one and two loop 
corrections liil' and the pole mass corrections. The theoretical value of m^ is still 
uncertain by about 3 GeV, and so the experimental bound of ruh > 114 GeV will 
(conservatively) be interpreted as a theoretical evaluation of my, > 111 GeV (and 
similarly ruh > 120 GeV as ruh > 117 GeV). (2) Large tan/3 NLO corrections to the 
b ^ S'j decay rate lirP are included. (3) Loop corrections to rrih and rrij- are included 



(which are important for large tan/5). (4) QCD renormalization group corrections in 
going below the scale Msusy are used. (5) In the relic density calculation, all slepton 
neutralino co-annihilation effects are included with an analysis valid for large tan (3. 
There are now several groups that have carried out this calculation li^' lij' li^, 
and we have checked that they are generally in good agreement. 

We have not assumed any specific GUT group relations, such as Yukawa 
unification or proton decay constraints (except for grand unification of the gauge 
coupling constants). Such relations depend on physics beyond the GUT scale about 
which little is known, and need not, for example, hold in string models, even when 
coupling constant unification is required. 

3 mSUGRA Model 

The mSUGRA model lif is the simplest possible model with universal soft breaking 
masses at Mq. It depends on four parameters and one sign: mo, the masses of the 
scalar particles at Mq] ^1/2, the gaugino masses at Mq; Aq, the cubic soft breaking 
mass at Mq; tan/3, the ratio of the Higgs VeVs, {H2) / {Hi) , at the electroweak scale; 
and the sign of fi, the Higgs mixing parameter in the superpotential (W^ = fiHiH2). 
We allow the above parameters to vary over the following ranges: < rriQ, mi/2 < 1 
TeV, 2 < tan/3 < 45, \Ao\ < 4mi/2. 



It has been known from the beginning B' M) that qSUGRA increases with 
tan p. The leading terms come from the chargino (xi ) diagrams where the muon 
can polarize into a chargino and a sneutrino. Expanding for {fi ± 7722)^ ^ M^ 



(which generally is the case for mSUGRA), one finds for this leading term liS 



,SUGRA 



a 



m. 



tan/9 



4:71 sin 9w 



m-±n / 1 _ !!!i 



M^^ 



1 + 3^ 



/i^ 



1-4'^ 



F(x) 



where m~± 



r?7,2 — 0.8mi/2, a^ = rnl/m^-± and F(a;) is a form factor arising from 



the loop integration. This implies that the sign of a' 



SUGRA 



is determined by the sign 



of /i VJ\) ^ and since experimentally, Eq. (2), one has that a^ has a positive anomaly, 
one has that 

/i > (9) 

This result has immediate consequences for dark matter detection. Thus 
the scattering of neutralinos by nuclei in dark matter detectors is governed by the 
cross section o'^o.p, and for /i < 0, cancellations can occur reducing these cross 

sections to below 10^^^ pb in large regions of the parameter space liS' Uj/. This is 
exhibited in Fig. 1. If /i were negative, one would require a detector well beyond 
anything that is currently planned, and so dark matter experiments would not be 
able to scan the full SUSY parameter space. However, with Eq. 9, one finds that 
cross sections generally are greater than 10"^'' liMP, accessible to future planned 
detectors such as GENIUS "* or Cryoarray cl^) . 
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Figure 1: cr^o_p for /i < 0, Aq = 1500 GeV, for tan/3 = 6 (short dash), tan/3 
(dotted), tan/? = 10 (solid), tan/? = 20 (dot-dash), tan/? = 25 (dashed) > . 



Eq. (9) also has important theoretical consequences for the mSUGRA 
model. Thus we first note that the combined experimental constraints on the Higgs 
mass and the b ^ s'-f branching ratio, implies that most of the allowed parameter 
space is in the co-annihilation region of the parameter space, i. e. mi/2 ~ 300 GeV. 
This means that rriQ is essentially determined in terms of mi/2 (for fixed Aq and 
tan/5) as can be seen in Fig. 2. In particular, this means that rriQ is an increasing 



function of mi/2. Now a 



SUGRA 



decreases with increasing mi/2 and increasing mo. 



Thus the lower bound of a^^*^^"^ of Eq. (7) then implies an upper bound on mi/2 
(since an increase in mi/2 cannot be compensated by a decrease of mo). Thus the 
combined effects of the data is to give now both lower bounds on 1711/2 (from rrih 
and b — > 37) and upper bounds on mi/2 (from a^), as illustrated in Fig. 3. If 
however, the a^ data had required /i < (instead of Eq. (9)) with the same lower 
bound on |a^^*^^^|, then the ruh and b -^ S'j constraints for this sign of /i would 
have essentially eliminated all the mSUGRA parameter space. Thus Eq. (7) also 
represents an experimental test of the validity of the mSUGRA model. 
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Figure 2: Corridors in the mo — m,i/2 plane allowed by the relic density constraint 
for tan/3 = 40, nih > HI GeV, /i > for Ao = 0, — 2?7ii/2,4mi/2 from bottom to 
top. The curves terminate at low mi/2 due to the b ^ S'y constraint except for the 
Aq = 4mi/2 which terminates due to the m^ constraint. The short lines through the 
allowed corridors represent the high mi/2 termination due to the lower bound on a^ 

ofEq. (7). 16) 



Since mn is an increasing function of rrii/2 and tan/5, the fact that rTii/2 
is bounded from above, and nih is bounded from below implies a lower bound on 
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Figure 3: mSUGRA contribution to a^ as a function of mi/2 for Aq = 0, n > 0, for 
tan/3 = 10, 30 and 40 (bottom to top). ^6) 



tan/3. We find at the 95% C. L. that for m/, > 114 GeV 

tan/5 > 7(5) for Aq = 0(-4mi/2) 
and for rrih > 120 GeV one would have 

tan/3 > 15(10) for Aq = 0(-4mi/2) 



(10) 



;ii) 



We see that the a^ anomaly favors large tan p. Fig. 3 also shows that the mSUGRA 
model cannot accommodate large values of a^^'"'^^, and if the final E821 results 
were significantly higher than 50 x 10~^°, this would be a signal for non-universal 
soft breaking. 

The strong constraints on mi/2 discussed above also sharpens considerably 
the predictions that mSUGRA has for the SUSY mass spectrum expected at accel- 
erators. To see what these are, we consider the 90% C. L. on the a^ anomaly, i.e. 
we require now a^^^^^ > 20 x 10~^° cl^ . In this case one finds that for Aq = 
that tan/3 > 10. Then for tan/3 < 40, mi/2 and mo are constrained as follows: 



mi/2 = (190 - 550) GeV; mo = (70 - 300) GeV 



(12) 



In Table 1 we give the range expected for the SUSY masses implied by 
these constraints. All the masses except the light stop, ti, are insensitive to Aq and 
the lower bound for m^^ can be lowered to 240 GeV by decreasing Aq to — 4mi/2- 



From this we see that the trilepton signal unfortunately would be out of reach of 
the Tevatron RUN II (tan/? and mi/2 are too large ED) and RUN II (with 30 fb"^) 
would only be able to see the light Higgs provided nih < 130 GeV cLy . A 500 GeV 
NLC would be able to detect the h, fi, and cover only part of the parameter space 
for the light selectron, ei. The LHC, however, would be able to see the full SUSY 
spectrum. 



Table 1. Allowed ranges for SUSY masses in GeV for mSUGRA assuming 
90% C. L. for a^ for Aq = 0. The lower value of mi^can he reduced to 240 GeV by 
changing Aq to -Ami/2. The other masses are not sensitive to Aq. ' ' 
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Figure 4: cr^o_p as a function of the neutralino mass m^o for tan/3 = 40, /i > for 
Aq = — 2mi/2, 4r?2i/2, from bottom to top. The curves terminate at small m^o due 
to the b —>■ s'-f constraint for Aq = and —2mi/2 and due to the Higgs mass bound 
(nth > 114 GeV) for Aq = Ami/2- The curves terminate at large m^o due to the 

lower bound on a^ of Eq. (7). ^ 

We saw above that the experimental requirement that a 



SUGRA 



be positive, 

i. e. that /x > 0, eliminated much of the parameter space that would not be accessible 
to dark matter detectors. In addition, the lower bound of Eq. (7) produces an upper 
bound on mi/2 and mo, and since the dark matter detection cross section, cr^o_p, 
decreases with increasing mo and mi/2, the lower bound on cr^o_p is raised for the 



remaining /i > part of the parameter space. Fig. 4 shows the expected cross 
sections for tan/3 = 40, /i > 0, rrih > 114 GeV, for Aq = — 2mi/2, 4mi/2, (bottom 
to top). We see that there is a significant dependence on Aq, and over the full range, 
cr^o_p > 6 X 10^^*^ pb. If the Higgs bound was raised to rrih > 120 GeV, then the 
lower bounds on r7ii/2 increase to 200 GeV, 215 GeV, and 246 GeV respectively, 
significantly further reducing the parameter space. 

If we reduce tan /? one would expect cr^o_p to fall. However, for lower tan/3, 
the upper bound on mi/2 becomes more constraining, eliminating more of the high 
^^1/25 ""^0 region, and thus compensating partially for the reduction of tan/3 . This 
is seen in Fig. 5 where cr^o_p is given for tan/3 = 10, Aq = (upper curve) and 
Aq = —4:1711/2 (lower curve). We have now that cr^o_p > 4 x 10~^° pb. Thus almost 
all the parameter space should now be accessible to future planned detectors such 
as GENIUS and Cryoarray. 
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Figure 5: cr^o_p as a function of m^o for tan/3 = 10, yU > 0, rrih > 114 GeV for 
Aq = (upper curve), Aq = —4mi/2 (lower curve). The termination at low m^o 
is due to the rrih bound for Aq = 0, and the b -^ S'-f bound for Aq = — 4r?2i/2. The 

termination at high m^o is due to the lower bound on a^ of Eq. (7). ^ 



4 Non-Universal Models 

We consider here the case where there is non-universal soft breaking masses at 
Mg for the Higgs bosons and the third generation of squarks and sleptons. This 
possibility produces interesting new effects. One may parameterize the masses in 



the following way: 



"^ii = "^0(1 + ^1)5 m^^ = ml{l + 62); 



m 



QL 



^l{^ + h)\ "^t« ="^o(l + '^4); m^=ml{l + 5^); 



Wa = "^0(1 + ^6); mil = ml{l + 5j). 



(13) 



where ttiq is the universal soft breaking mass for the first two generations, qi 
(iL,bL), II = {t^l.t'l), and we assume that —l<5i< +1. 
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Figure 6: Effect of a non-universal Higgs soft breaking mass enhancing the Z^ s- 
channel pole contribution in the early universe annihilation, for the case of 62 =1, 
tan/5 = 40, y4o = wii/2, /i > 0. The lower band is the usual fi coannihilation region. 
The upper band is an additional region satisfying the relic density constraint arising 
from increased annihilation via the Z^ pole due to the decrease in /i^ increasing the 
higgsino content of the neutralino. > 

While this model has a number of additional parameters, one can under- 
stand the physics of what is implied from the following considerations. Electroweak 
symmetry breaking is governed by the ix parameter at the electroweak scale, and for 
low and intermediate tan/5 one has 
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/^ 



t2-l 
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Dn 
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0'). + -7^ 



ml + universal parts + loop corrections. 



(14) 



Here t = tan/?, and Dq = 1 — {mt/200 GeVsin/?)^ = 0.25 . {Dq arises from running 
the RGE from Mq to the electroweak scale.) One sees that for tan/3 > 3, the 



universal part of /i^ depending on ml is quite small, and so this term is sensitive 
to the amount of non- universal soft breaking that might be present. Now much 
of the physics is governed by fi^. If /x^ is decreased, then the higgsino part of Xi 
will increase. This has several effects. First it increases the x\ ^ x\~ ^ coupling 
which allows then the opening of a new neutralino annihilation channel through 
an s-channel Z pole, in the relic density calculation. This allows a new region of 
allowed relic density at high tan/5 and high mo. Second, since cr^o_p depends on the 
interference term between the higgsino and gaugino parts of the neutralino, the dark 
matter detection cross section will be increased with increasing higgsino content. To 
illustrate these effects, we consider two cases 
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Figure 7: cr^^o-p o-s a function of mi/2 (m.^o = 0.41711/2) for tan (3 = 40, fi > 0, 
TJifi > 114 GeV, Aq = mi/2 for 62 = I- The lower curve is for the fi — Xi co- 
annihilation channel, and the dashed band is for the Z s-channel annihilation allowed 
by non-universal soft breaking. The curves terminate at low mi/2 due to the b ^ S'y 
constraint. The vertical lines show the termination at high mi/2 due to the lower 

bound on a^ of Eq. (7). ^ 

(1) 62 = 1, 6, = 0, ^ ^ 2. 

Here one sees that the tUq term contributes negatively to /i^. The effect 
of this is shown in Fig. 6 where the region allowed by the relic density constraint 
in the rrio — mi/2 plane for tan/5 = 40 , Aq = mi/2 is given. The lower corridor, 
coming from the stau-neutralino co-annihilation is much the same as the Aq = 
corridor in Fig. 2. The upper channel, at higher niQ, is due to the s-channel Z- pole 
annihilation. Fig. 7 shows the neutralino proton cross sections for corresponding 
allowed regions as a function of mi/2 {m^o = 0.4mi/2), for tan/5 = 40. In spite of 
the fact that niQ is quite high for the Z s-channel corridor, the cross sections is quite 



large, and at least part of the this region should be accessible to GENIUS -TF, and 
CDMS in the Soudan mine. This increase is due to the fact that the non-universality 
has lowered ii^. 
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Figure 8: Allowed regions in the mQ — mi/2 plane for the case tan/3 = 40, Aq = mi/2, 
fi > 0. The bottom curve is the mSUGRA fi coannihilation band of Fig. 2 (shown for 
reference). The middle band is the actual fi coannihilation band when 6iq = —0.7. 
The top band is an additional allowed region due to the enhancement of the Z^ s- 
channel annihilation arising from the nonuniversality lowering the value of //^ and 
hence raising the higgsino content of the neutralino. For mi/2 ~ 500 GeV, the two 

bands overlap. ) 



(2) (5io(= 5s = S^ = 5,) = -0.7. 

This is a model that might arise in a SU{5) GUT, where only the particles 
of the 10 representation in the third generation have non-universal soft breaking 
masses. Again, the non- universalities lower the value of /i^, allowing for a Z s- 
channel corridor in the relic density analysis. In addition, because the soft breaking 
mass of the stau is reduced, the co-annihilation channel occurs at a higher value of 
mQ. This is shown in Fig. 8, for tan/3 = 40, where the bottom band, the usual 
mSUGRA CO- annihilation channel, is now moved up to much higher nio. (This co- 
annihilation and the Z-channel annihilation corridors merge for mi/2 ~ 500 GeV.) 
The corresponding cross sections are shown in Fig. 9. Again part of this parameter 
space should be accessible to the next round of dark matter detectors. 
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Figure 9: cr^o_p as a function of mi/2 for tan/3 = 40, fi > 0, Aq = mi/2 and 
rrih > 114 GeV. The lower curve is for the bottom of the fi — x^ co- annihilation 
corridor, and the upper curve is for the top of the Z channel hand. The termination 
at low 7Tii/2 is due to the b ^ S'-f constraint, and the vertical lines are the upper 

bound on mi/2 due to the lower bound of a^ of Eq. (7). ) 

5 Conclusions 

We have examined here the 2.6 a deviation of a^ = (gf^— 2)/2 from its predicted Stan- 
dard Model result under the assumption that the effect is real and can be understood 
in terms supergravity GUT models with R-parity invariance. For mSUGRA models, 
the combined constraints from the a^, nih, and the b —>■ sj decay constraints and 
the relic density of Xi dark matter greatly limit the SUSY parameter space, greatly 
tightening the predictions the theory can make. The lower bound on a^ produces 
an upper bound on mi/2, and then the m^ and b ^ s'y constraints produce lower 
bounds on mi/2 and tan/3. We find that for rrih > 114 GeV, that tan/3 > 7(5) for 
Aq = 0(-4mi/2) and for rrih > 120 GeV, tan/3 > 15(10) for Aq = 0(-4mi/2). The 
lower bound on r7ii/2 pushes most of the parameter space into the co-annihilation do- 
main effectively determining mo in terms of mi/2- However, there is still a significant 
dependence on Aq and tan/3. 

The restrictions on the allowed range of mi/2 and ttiq allow one now to 
make significant predictions concerning accelerator reaches for the SUSY particles. 
Thus we find at 90% C. L. bounds on a^ that the Tevatron Run II should be able 
to see the light Higgs (if nih < 130 GeV), but no other parts of the SUSY spectrum. 
A 500 GeV NLC should be able to see the h and fi, and scan part of the parameter 
space where the ei is expected. The LHC should be able to detect the full SUSY 



spectrum. Further, future planned dark matter detectors, GENIUS and Cryoarray, 
should be able to sample almost all of the SUSY parameter space. For other analyses 



within the mSUGRA framework see ^^K 

Non-universal SUGRA models can allow new regions in parameter space 
to open due to possible annihilation in the early universe through s-channel Z poles. 
For these types of non- universalities, which lower the value of fi'^, the neutralino 
proton cross sections become larger, and part of these effects can possibly be tested 
with current detectors. Since non-universalities are phenomena determined by post- 
GUT physics, the dark matter detectors could allow one to learn about physics 
beyond the GUT scale. 

The current anomaly in a^ is a 2.6 a effect. However, the additional data 
that now exists, and the expected reduction of the theoretical errors should allow a 
determination in the near future of whether the effect is real. 
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